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We employ one-particle-exchange method to study D — > PI/ decays in D ^ Kp, 
C^ij 7r.fr*, irp processes. Taking into account a strong phase and considering nonfactorizable 

effect, we can get good results consistent with the experimental data. Nonfactorizable 
effect is not always large, but in some cases, the nonfactorizable effect is necessary 
to accommodate the experimental data. Strong phase is approximately SU(3) flavor 
^ ' symmetric. 



Abstract 



1 Introduction 

To understand the quark mixing sector of the standard model(SM) and search for new physics 
beyond the SM, one needs to study the decays of heavy mesons and calculate precisely the 
transition matrix elements of the heavy mesons decays. The short distance effects due to 
hard gluon exchange can be calculated reliably and the effective hamiltonian and factoriza- 
tion approach has been constructed [[I], 0], thus a lot of results which fit the experimental 
data well have been obtained by factorization approach. However, there are still many decay 
modes which can not be accommodated by factorization approach. In fact, the quarks in 
heavy mesons are bound by strong interaction which is described by nonperturbative QCD. 
* Supported in part by National Natural Science Foundation of China. 
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After weak decays of heavy mesons, the final particles can rescatter into other particle states 
through nonperturbative strong interaction || ||, this is called final state interaction (FSI). 
Many authors have studied FSI effects and found that FSI effects may play a crucial 
role in some decay modes [[5], ||. Therefore it is necessary to study heavy meson two-body 
weak decays beyond the factorization approach. Since the FSI process refers to the soft 
rescattering process which is controlled by nonperturbative QCD and can not be reliably 
evaluated with well-established theoretical frame, we have to rely on phenomenological mod- 
els to analyze the FSI effects in certain processes. One can model this rescattering effect as 
one-particle exchange process |7|, |8|]. There are also other ways to treat the nonperturbative 
and FSI effects in D decays, the readers can refer to Ref.||. In this paper, we study some 
channels of D — > PV decays. We use the one-particle-exchange method to study the final 
state interactions in these decays. The magnitudes of hadronic couplings needed here are 
extracted from experimental data on the measured branching fractions of resonance decays. 
In addition, we consider a strong phase for the hadronic coupling (JToJ which is important 



for obtaining the correct branching ratios of D — > PV decays. We also take into account 



some possible nonfactorizable effect flll], |12| , which is needed for some decay mode from 
the phenomenological point of view. The coupling constants extracted from experimental 
data are small for s-channel contribution and large for t-channel contribution. Therefore the 
s-channel contribution is numerically negligible in D — > PV decays. We can safely drop the 
s-channel contribution in this paper. 

The paper is organized as follows. Section II presents the calculation in naive factorization 
approach. Section III gives the main scheme of one-particle-exchange method. Section IV is 
devoted to the numerical calculation and discussions. Finally a brief summary is given. 

2 Calculations in the factorization approach 

The low energy effective Hamiltonian for charm decays is given by fl3|| 

G 

Ti-eff = -^{V us V* s [C 1 (sc) v _ A (us)v_ A + C 2 (ss) v _ A (uc)v_ A ] 
+VudV* d [C 1 (dc) v ^ A (ud) v ^ A + C 2 {dd) v _ A (uc) v „ A ] 
+V ud V* s [C 1 (sc) v _ A (ud) v _ A + C 2 {sd) v ^ A {uc) v _ A }} + h.c. , (1) 

where C\ and C 2 are the Wilson coefficients at m c scale. We need not consider the con- 
tributions of the QCD and electroweak penguin operators in the decays of D — > PV, since 
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their effects are small in D decays. The values of C\ and C 2 at m c scale are taken to be 
C l = 1.216, C 2 = -0.415 H 

In the naive factorization approach, the decay amplitude can be generally factorized into 
a product of two current matrix elements and can be obtained from eq.(]l|) 
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(2) 



where the parameters a\ and a 2 are taken as [0 

1 



ax = cx(fi) + c 2 (/i) ( — + x(m) 



a 2 



c 2 (m) + ci(/x) ( + 



(3) 



with the color number N c = 3, and is the phenomenological parameter introduced for 
taking care of nonfactorizable effects. The parameters in calculation are: 1) the form factors, 
Ff n (0) = 0.69, F? K (0) = 0.76, Aq P (0) = 0.67, Ag K *(0) = 0.73 §; 2) the decay constants, 
U = 0.133GeV, f K = 0.158GeV, f p = 0.2GeV, and f K . = 0.221GeV. 

For q 2 dependence of the form factors, we take the BSW model @, i.e., the monopole 
dominance assumption: 

*i(0) „ M0) 



Fx(q 2 



A (q 2 



1 — g 2 /m^_' ' ' \ — q 2 /m 2 ] -'' 

where mi-,m - is the relevant pole mass. 

The decay width of a D meson at rest decaying into PV is 



(4) 



T(D -> PV) = — \A{D -> PV)\- 
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D 



(5) 



where the momentum of the final state particle is given by 

[(m 2 D - (mi + m 2 ) 2 )){{m 2 D - (m x - m 2 ) 2 )] 1/2 



2m 



D 



(6) 



where m\,mi are the masses of final state particles. The corresponding branching ratio is 

T(D -> PV) 



Br(D -> Py) 



(7) 



tot 



Table 1: The branching ratios of D — > PI/ obtained in the naive factorization approach and 
compared with the experimental results. 



Decay mode 


Br (Theory) 

(x = o) 


Br (Theory) 
(ai = 1.26,02 = -0.51) 


Br (Experiment) 


D° -> K°p° 


3.92 x lO" 3 


5.92 x 10~ 3 


(1.21 ±0.17) x 10~ 2 




10.66 x 10~ 2 


11.45 x 10" 2 


(10.8 ±0.9) x 10~ 2 


p+ -> if°p + 


17.35 x 10~ 2 


16.91 x 10~ 2 


(6.6 ±2.5) x lO" 2 


L>° -> 7T°K$ 


1.37 x 1(T 2 


2.08 x 10~ 2 


(3.1 ±0.4) x 10~ 2 


D° -> 7T + ^- 


3.06 x 10~ 2 


3.29 x 10~ 2 


(5.0 ±0.4) x 10~ 2 


£)+ _> tt+K* 


8.72 x lO- 3 


3.66 x 10~ 2 


(1.90 ±0.19) x 10- 2 


P+ -> 7T + p° 


8.12 x 10~ 3 


9.33 x 10~ 3 


(1.05 ±0.31) x 10~ 3 


P U -> TT+p- 


1.36 x 10~ 3 


1.46 x 10~ 3 




P° -> 7T p° 


7.17 x 10" 4 


1.09 x 10~ 3 




L>° -> 7T~p+ 


4.48 x 10" 3 


5.89 x 10" 3 




P+ -> 7T U p + 


1.78 x 10~ 2 


1.98 x 10~ 2 





The numerical results of the branch ratios of D decays are given in Table 1. The second 
column is for the case x(/x) = which means there is no nonfactorizable contribution. When 
x(/i) = 0, the parameters ai = 1.216 and a 2 = —0.415. The parameters in the third column 
a\ = 1.26, ci2 = —0.51 are phenomenologically used in many references [14], which is relevant 
to taking non-zero parameter xit 1 )- 

Comparing the results of the naive factorization in the second and third column of Table 
1 with the experimental data, one can notice that, even considering some nonfactorizable 
contribution, some of the results from the naive factorization approach deviate significantly 
from the experimental data. 

3 The one particle exchange method for FSI 

From Table 1, we can see that the calculation from naive factorization approach is in dis- 
agreement with the experimental results for the branching ratios of D — > PV decays. The 
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reason is that the physical picture of naive factorization is too simple, in which nonper- 
turbative strong interaction is restricted in single hadrons, or between the initial and final 
hadrons which share the same spectator quark. If the mass of the initial particle is large, 
such as the case of B meson decay, the effect of nonperturbative strong interaction between 
the final hadrons most probably is small because the momentum transfer is large. However, 
in the case of D meson, its mass is not so large. The energy scale of D decays is not very 
high. Nonperturbative effect may give large contribution. According to the one-particle 
exchange method, there are s-channel and t-channel contribution to the final state interac- 
tion ]7|, ||. The diagrams of these nonperturbative rescattering effects can be depicted in 
Figs.|l| and 0. The first part D — > P1V2 or D — > V1P2 represents the direct decay where the 
decay amplitudes can be obtained by using naive factorization method. The second part 
represents rescattering process where the effective hadronic couplings are needed in numer- 
ical calculation, which can be extracted from experimental data on the relevant resonance 
decays. 




Figure 1: s-channel contributions to final-state interaction in D — > PV decays. 

There are many resonances near the mass scale of D meson, it is possible that nonper- 
turbative interaction is propagated by these resonance states, such as, if* (892), if* (1430), 
/o(1710), if*(1680), if*(1020), 0(1680), vr(1300), etc. For s-channel the correct quantum 
number of the resonance should be J p = 0~ (in charged D decays) . In neutral decay modes, 
the reasonance should be with quantum number J PC = l ". For D° — > ir°p , only 7r(1300) 
has the correct quantum number f!5 |. Fig.[l| is the s-channel contribution to the final state 



interaction in D° — > 7r°p°. Here V\ and P2 are the intermediate mesons. Because the cou- 
pling of 7r(1300) with 7r°p° is too small [|Kj, we can ignore the s-channel contribution in the 
numerical analysis. Fig.^| shows the t-channel contribution to the final state interaction. Pi, 
V 2 and Vx, P2 are the intermediate states from direct weak decays. They rescatter into the 
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Figure 2: t-channel contributions to final-state interaction in D 
change, (a) D -► PiV 2 -» PV, (b) L> -> FiP 2 -> PV. 



due to one particle ex- 



final states by exchanging one resonance state P. In this paper the intermediate states are 
treated to be on their mass shell, because their off-shell contribution can be attributed to the 
quark level. We assume the on-shell contribution dominates in the final state interaction. 
The exchanged resonances are treated as a virtual particle. Their propagators are taken as 
Breit-Wigner form, 

72 2~T~ — F — ' ( 8 ) 

Ar — m z + iml t t 

where T to t is the total decay width of the exchanged resonance. 

We consider the t channel contribution. For the t-channel contribution, the concerned 
effective vertex in Fig.[5] is VPP, which can be related to the V decay amplitude. Explicitly 
the amplitude of V — ► PP can be written as 



Typp = 9vpp e • (pi - P2), 



(9) 



where pi and P2 are the four-momentum of the two pseudoscalars, respectively. To extract 
gvpp from experiment, one should square eq.(Bf) to get the decay widths, 



Yiy^PP) = ii-i^ppf 

6 87T 



m v — 2m l — 2m 2 + 



2\2 



V 



rn 



2 • 



(10) 



where mi and m 2 are the masses of the two final particles PP, respectively, and | p | is the 
momentum of one of the final particle P in the rest frame of V. From the above equations, 
one can see that only the magnitudes of the effective couplings | gypp | can be extracted from 
experiment. On the quark level, the effective vertex should be controlled by nonperturbative 
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Figure 3: The effective coupling vertex on the hadronic level 

QCD. It is reasonable that a strong phase can appear in the effective coupling, which is 
contributed by strong interaction. Therefore we can take a strong phase for each hadronic 
effective coupling |TIJ. In the following, the symbol g will only be used to represent the 



magnitude of the relevant effective coupling. The total one should be ge l9 , where 6 is the 
strong phase. For example, the effective couplings will be written in the form of gyppe %6vpp ■ 
The t-channel contribution in Fig.^(a) is 



Pl ' V2 2 J (2ir) 3 2E 1 J (2vr) 3 2P 2 

I e i(Oi+0 2 ) 

where F(k 2 ) = (A 2 — m 2 ) / (A 2 — k 2 ) is the form factor which is introduced to compensate the 



off-shell effect of the exchanged particle at the vertices [16[ . We choose the lightest resonance 
state as the exchanged particle that gives the largest contribution to the decay amplitude. 
After a few steps of integration to the above equation, we get 

FSI /-i djcosO) »e'(ft+fc> 2 2 

A Pi,V2 ~ / ~7y Pi A-i fi'iyi 2 , • -p b \ k ) 92 111 , (12) 

where 



H x = m 2 ft A [-(E 1 E 4 + |pi||p 4 | cos6>) 

+ -^(M 2 - m\ - m 2 2 )(E 2 E 4 - |p 2 ||p 4 | costf)] 

x[— (\n\Et - E 3 \n\)}, (13) 
nT>3 

and Xi represents the relevant direct decay amplitude of D decaying to the intermediate 
pair P x and V 2 divided by (P^V - A)J0) (V 2 \(V - Ay\D), 

A(D -> P X F 2 ) 



(-PiK^ - A)J0)(V 2 \{V - Ay\D) ' 
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The t-channel contribution in Fig.^(b) is 



% e 



i(0i+e 2 ) 



*9i e x ■ (p 3 - k)— s— -= F(k ) g 2 e 4 ■ (p 2 - fc), (14) 



A; 2 — m 2 + imr 



to; 



and we obtain 



/•! d(cos0)_ . i e^+o^ 



where 



#2 = mihFi[-M D Ez + ExM D (ExE 3 - costf)] 

x— (|p 4 |£ 2 -£ 4 |p 2 |cos#), (16) 
m 4 

and X 2 represents the relevant direct decay amplitude of 2? decaying to the intermediate 
pair Vi and P 2 divided by - A) ^0} (P 2 \(V - A)*\D), 

x A(D -> V^) 

2 " (K|(V-AV|0)<^|(V-^|D>- 

4 Numerical calculation and discussions 

To calculate FSI contribution of D decays with the eq.(|T2D and eq. (|15D , we need to know 
which channels can rescatter into the final states. For D — > Kp, itK*, Tip, from Figs.f|~^], one 
can see that D -> ixK* ^ Kp, D ^ Kp ^ Kp, D -> vr^* -> ttA"*, D -> pif -> tt/^*, L> -> 
7rp — > 7rp, D — > KK* — > 7rp can give the largest contributions, because these intermediate 
states have the largest couplings with the final states and the masses of exchanged meson 
are small which give the largest t-channel contributions. When calculating the contribution 
of each diagram in Figs.f|~|5], we should, at first, consider all the possible iso-spin structure 
for each diagram in Figs.|]~§], and draw all the possible sub-diagrams on the quark level. 
Second, write down the iso-spin factor for each sub-diagram. For example, the uu component 
in one final meson p° contributes an isospin factor and the dd component contributes 
— A=. For the intermediate state p°, the factor and — 4j should be dropped [|IJ]. Third, 



sum the factors of all the possible sub-diagrams of each diagram to get the iso-spin factor 
for each diagram on the hadronic level. For example, in the diagram (a) of D° — > K°p°, the 



S 




Figure 4: Intermediate states in rescattering process for D — » Kp decays. 

iso-spin factor of one sub-diagram is ^= and is — ^ in another sub-diagram, so the factor of 
the diagram (a) of D° —>■ K°p° is zero. 

From eq. fli~2l) , eq.(|l5|) and considering Figs.|]-^|], we can calculate the amplitudes of 
D — > PV decays. In this paper we consider D — > Kp, ttK* and up decays. There should be 
some input parameters in our calculation, such as, the transition form factors for D decays, 
decay constants of the final mesons, the phenomenological nonfactorizable parameter x(p)i 
the off-shell compensating parameter A in function F(k 2 ) introduced in eq. QTTl) , the effec- 
tive couplings of relevant hadronic states and the relevant strong phases for these effective 
couplings. For the transition form factors and decay constants, we take 1) the form factors, 
F(>"(0) = 0.69, Ff K (0) = 0.76, Aq P (0) = 0.67, A° K *(0) = 0.73 @; 2) the decay constants, 
U = 0.133GW, f K = 0.158GW, f p = 0.2GeV, and f K * = 0.221GeV. We should be careful 
for these parameters, because except for the decay constants f n and fx, etc., the values of 
the transition form factors have not been known exactly yet. We have to take them from 
model-dependent calculations. For the phenomenological nonfactorizable parameter x{p)i 
at first, we tried to proceed by taking x{p) = 0> which means that nonfactorizable con- 
tribution is neglected. We find that if nonfactorizable contribution is neglected, no matter 
how the other parameters ( the strong phases and A ) are tuned, we can not reproduce the 
experimental data for all the D — > PV decays simultaneously. So we have to keep it as 
an phenomenological parameter which will be determined later. The hadronic effective cou- 
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Figure 5: Intermediate states in rescattering process for D — * ttK decays. 

plings involved in this study are g p7Tn and gK*Kn, which can be determined from the center 
value of the measured decay width of p — > ixix and K* — > Kit WB. We obtain g pn7T = 6.0, 
9k* Kit = 4.6. The parameter A in the off-shellness compensating function F(k 2 ) introduced 
in eq. fllTD is not an universal parameter. It is pro cess- dependent in general. However, in this 
paper we use one value for A in all the possible channels of D — > PV decays. We assume A 
is in the range from 0.5GeV to l.OGeV, which is the range of the masses of the final state 
particles, p and K*, etc.. We scanned all the possible value for x(f^) an d A, and find that 
if we take x(aO = 0-16 and A = 0.7GeV, we can reproduce the experimental data of all the 
detected D — > PV decay modes well. x(/-0 = 0-16 means the nonfactorizable contribution is 
not large. A = 0.7GeV is in the mass range of the final state particles. In the following we 
give the decay amplitudes of some D — > PV decay modes as function of the strong phases 
0k*Ktt, Qp-w-K and 9 p kk by taking xip) = 0-16 and A = 0.7GeV, 

A(D° -> K°p°) = -6.768 x 10~ 7 + 1.40508 x 10- 6 ie i28 » KK , 

A(D° -> K-p + ) = 4.09 x 10" 6 - 2.2512 x 10~ 7 ie i2 ^ K - 6.739 x io- 7 ie < P**** +fl ' n " t ), 

A(D + -> K°p + ) = 3.2579 x 10~ 6 + 1.14838 x 10- 6 ie i26 » KK + 6.114 x 10~ 7 'ie i(0Jf * , 

A(£> -> 7r°i? *) = -1.1837 x 10" 6 + 1.128 x 10 _6 ze i2 ^*^, 

-> 7r + lT~) = 2.352 x 10~ 6 - 5.501 x 10~ 7 ie i2 ^*** - 4.51 x lO^ie^*^ 9 "^, 

A(D + -> 7T+i? *) = 1.2056 x 10" 6 + 3.7368 x lO^ie* 2 **"'** + 1.81603 x i(r 6 ie <(ejc * Jt * + *'"" r >, 
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Figure 6: 
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Figure 6: Intermediate states in rescattering process for D — > up decays. 

-6.728887 x 10~ 7 - 9.61 x 10- 7 ie i26 » KK + 2.2334286 x 10" 7 i e i{eK * K - +e " KK \ 
-4.8216 x HT 7 + 1.137 x \^ i e ^K*K^ pKK ) +2 .\x l(T 7 ie i2f ^ , 
-2.608 x 10~ 7 + 1.821 x \^ ie i(dK * K ^^ KK \ 
-8.736 x 10~ 7 + 6.879 x iO- 8 ie i( - e ^ K - +e " KK) + 2.08 x 10~ 7 ie i2e "^ , 
-1.0597 x 10~ 6 - 9.6106 x l^ie 126 "^ + 2.21 x nr^e^*^^. (17) 

The phases of the effective hadronic couplings 9k* Km 9 p7Tn and Q p kk can not be known 
from direct experimental measurement or from any nonperturbative calculations because 
there are no any such kind of calculations yet. We only know that the values of 0K*Kn, 
6 pwn and 9 p kk should not differ too much according to SU(3) flavor symmetry. We tried 
some values for these phase parameters, and find that the ranges which can reproduce the 
experimental data of the measured D — > PV decays are not very narrow. To show the 
situation that the experimental data are accommodated, we give the numerical results for 
0k*k* = 51.0°, 9 p kk = 51.0° and 9 pn7T = 57.3° in table 2. 

Table 2 shows that the contribution of FSI is strongly channel dependent. For example, 
For D° — > K°p°, the braching ratio in naive factorization is 3.922 x 10~ 3 , while the braching 
ratio including FSI is 1.25 x 10~ 2 . we can see that FSI contribution in D° — > K°p° is large, 
but FSI contribution in D° — > K~p + is small. The reason for the difference is that the 
external rescattering diagrams for D° — > K°p° and D° — > K~ p + are different. Without the 
contribution of FSI, predictions of naive factorization for most detected D — > PV decays 
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Table 2: The branching ratios of D — > PV. 



Decay mode 


.Factorization 


ractonzation + rbl 


Experiment 


L> — > K u p 


o no - , i a — 3 

3.92 x 10 


1 or - , i n — 2 

1.25 x 10 


( 1 O 1 I A 1 r 7\ - , 1 A — 2 

(1.21 ± 0.17) x 10 


nO r^— + 
L> — > ii 


1 n or* - , 1 a — 2 
10. DO X 10 


11 1 ^ , 1 n — 2 

11.1 x 10 


/ 1 A O I AA\ - , 1 A — 2 

(10.8 ± 0.9) x 10 


iJ T — > K v p^ 


1 ( .oo x 10 


/.01 x 10 


l c c _i_ o K \ x^ 1 n — 2 
(D.D ± 2.5) X 10 


D° -> 7T°K * 


1.37 x 10~ 2 


2.72 x 10~ 2 


(3.1 ±0.4) x 10- 2 




3.06 x 10~ 2 


5.22 x 10~ 2 


(5.0 ±0.4) x 10~ 2 




8.72 x 10~ 3 


1.93 x 10~ 2 


(1.90 ±0.19) x 10~ 2 


-> 7T+p U 


8.12 x 10" 3 


1.3 x 10" 3 


(1.05 ±0.31) x 10~ 3 


L>° -> 7T + p- 


1.36 x lO" 3 


3.6 x 10~ 3 




D U -> 7tV° 


7.17 x 10~ 4 


1.1 x 10" 3 




D U -> 7T-p+ 


4.48 x 10~ 3 


7.3 x 10- 3 




-> 7T°p + 


1.78 x 10~ 2 


3.1 x 10~ 3 





are seriously in disagreement with the experimental results. After including FSI, the results 
can accommodate the experimental data well. For the other decay modes D° — > n + p~, vr°p°, 
7r~p + and D + — > 7r°p + , their branching ratios have not been detected in experiment yet. In 
our model, they are all predicted to be at the order of O(l0~ 3 ). For D° — > n + p~ and n~p + , 
the effect of FSI is enhancement. While for D° — > 7r°p° and D + — > 7r°p + , suppresses 
the prediction of naive factorization. 

Before the end of this section, some comments should be given: there are many uncer- 
tainties in the input parameters which may change the above result numerically, such as, 
the D decay transition form factors and some decay constants which have not been known 
exactly yet, they need to be measured from leptonic and semileptonic decays of D mesons 
which are quite possible in CLEO-C program in the near future. The other sources which 
may cause uncertainties are the shape of the off-shell compensating function F(k 2 ), or in 
more general the effective hadronic couplings in the off-shell region, the strong phases of the 
effective couplings, and the nonfactorization parameter x(/-0> both of which are needed to 
be studied in some nonperturbative methods based on QCD in the future. 

5 Summary 

We have studied some channnels of D — > PV decays. The total decay amplitude includes 
direct weak decays and final state rescattering effects. The direct weak decays are calculated 
in factorization approach, and the final state interaction effects are studied in one-particle- 
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exchange method. The prediction of naive factorization is far from the experimental data. 
After including the contribution of final state interaction, the theoretical prediction can ac- 
commodate the experimental data well. The strong phases of the effective hadronic couplings 
are neccessary to reproduce experimental data. 
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